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Abstract 

Detecting a stationary, stochastic gravitational wave signal is complicated 
by impossibility of observing the detector noise independently of the signal. 
One consequence is that we require at least two detectors to observe the sig- 
nal, which will be apparent in the cross-correlation of the detector outputs. A 
corollary is that there remains a systematic error, associated with the possible 
presence of correlated instrumental noise, in any observation aimed at esti- 
mating or limiting a stochastic gravitational wave signal. Here we describe 
a method of identifying this systematic error by varying the orientation of 
one of the detectors, leading to separate and independent modulations of the 
signal and noise contribution to the cross-correlation. Our method can be 
applied to measurements of a stochastic gravitational wave background by 
the ALLEGRO/LIGO Livingston Observatory detector pair. We explore — 
in the context of this detector pair — how this new measurement technique is 
insensitive to a cross-correlated detector noise component that can confound 
a conventional measurement. 

PACS numbers: 04.80.Nn, 04.80.-y, 95.55.Ym 
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I. INTRODUCTION 



While not yet detected, we can be sure that the Earth is bathed in a stochastic gravi- 
tational wave background. Contributions to that background arise from the confusion limit 
of large numbers of distant, conventional gravitational wave sources (e.g., binary systems), 
early universe physics (e.g., parametric amplification of fluctuations during an inflationary 
phase), and possibly other "exotic" physics (e.g., cosmic strings) Jll. One of the goals of the 
ground-based gravitational wave detectors now operating or under construction is to 
detect or place limits on the amplitude and spectrum of this stochastic gravitational wave 
background. 

A single, isolated gravitational wave detector cannot distinguish between instrumental 
noise and a weak, stationary cosmic gravitational wave background radiation. At least two 
detectors are needed, in which case the stochastic signal will be apparent in their cross- 
correlation. More particularly, the two detectors must i) have an overlapping frequency 
response, ii) have a separation shorter than the wavelength of their overlapping response, 
and in) both sample the same polarization state of the incident radiation. The LIGO Liv- 
ingston Observatory (LLO) interferometric detector |p!0|,|6| and the Louisiana State University 
ALLEGRO cryogenic acoustic detector P JTT|JT^ , separated by 42.3 km, constitute such a 
detector pair, capable of providing an experimental bound on the stochastic gravitational 
wave background at approximately 900 Hz. 

Unfortunately, a weak, stationary stochastic signal cannot be distinguished from a sim- 
ilarly weak, stationary noise background that is correlated between the two detectors. In- 
strumental noise arising from the environment may lead to such a correlated detector noise; 
consequently, a convincing case must be made that no terrestrial noise source is responsible 
for any observed correlation. This is a daunting experimental challenge. Any technique that 
can improve our ability to discriminate between a stochastic gravitational wave signal and 
a weak correlation of terrestrial origin should be pursued. 

Since the signal contribution to the cross-correlation depends on the relative orientation 
of the detectors (which determines their sensitivity to the two different polarization states), 
changing the orientation of one of the detectors will modulate the signal contribution to 
the cross-correlation in a predictable way, allowing us to distinguish the signal from the 
correlated noise and leading to a significantly improved estimate or bound on the in-band 
amplitude of a stochastic gravitational wave background. The technique of introducing a 
controlled signal modulation in order to identify and eliminate systematic environmental 
effects was conceived by Dicke |13| as the switching radiometer during the development of 
radar. 

The ALLEGRO group has used a planned relocation of their cryogenic detector to a 
new laboratory in order to implement the capability to re-orient the detector between data- 
taking periods This capability now allows for a modulation of the gravitational wave 
contribution to the detector noise cross-correlation in exactly the manner described. Here 
we discuss the details of the modulation and how it can be used to improve the reliability of 
the estimate or limit that we can place on the amplitude of a stochastic gravitational wave 
background near 900 Hz. 

The basic scheme for detecting a stochastic gravitational wave background using two or 



more gravitational wave detectors was described in |15| and further elucidated in |16-18 
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In Section |T| we review that work and extend it to the case where the two detectors share a 



cross-correlated noise component. In Section |T| we apply these results to the specific case of 



the ALLEGRO/LLO detector pair, providing the relevant geodetic and physical parameters 
that are needed to quantify the magnitude of the correlation and describing how to extract 
the correlation from observations. In Section |V| we present results of numerical calculations 
based on the initial LIGO instrumentation and its proposed upgrade. 

II. THE CROSS-CORRELATED DETECTOR OUTPUT 

A. The cross-correlation statistic 

The output of each detector — ALLEGRO or LLO — is a single time series, which is the 
sum of instrumental noise and a projection of the incident gravitational wave strain. Denote 
the output of LLO as sl and the output of ALLEGRO as sa and define the correlation of 
sl and sa over an integration time Tint by |T5Hl^ 



C{^Ia,^l) ■■= {sa,Sl) (la) 

/Tint/2 /"Tint/S 
dt / dt' SAit)sL{t')Qit - t';nA,nL), (lb) 
-Tint/2 J-Tint/2 



where 



Qa '■= Angles describing the orientation of ALLEGRO (Ic) 
fli '■= Angles describing the orientation of LLO. (Id) 

We discuss the choice of integration kernel Q (r; ft a, ^l) below. 

To evaluate the expectation value of C in the presence of signal and noise, write 

s* := + K, (2) 

where is the noise and h^^ is the signal in detector *. The signal is, of course, statistically 
independent of the noise; consequently, 

/i := C (3a) 

:= {Ha, hi) + {riA, til) , (3b) 

where the overbar represents an ensemble average. Introducing also the variance of C, 

a^:=C^-C\ (3c) 
we define the dimensionless signal-to-noise ratio (SNR) 

P ■= Pl^y- (4) 
As described above, the integration kernel in equation TO is at our disposal. If the detector 



and signal noise spectral densities are known, then Q can be chosen to maximize the signal 
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to noise ratio. Previou s work |1T6| - |T8| has always assumed that the contribution of the noise 
cross-correlation (e.g., {nA,nL)) to the mean /i vanishes. In Appendix ^ we consider the 
more general case where the noise contribution to the ensemble mean cross-correlation is 
non-zero. In this case the kernel Q that maximizes the SNR can be conveniently expressed 
in the frequency domain as 



where 



P(SA{f)SL{f) + SAL{fy 



(5a) 



Sa 
Sl 
Sal 

GW,0 



(ALLEGRO noise power spectral density) 

(LLO noise power spectral density) 

(ALLEGRO/LLO noise cross spectral density) 

/ expected stochastic signal spectrum expressed 

I fraction of the closure density in logarithmic frequency 



(5b) 
(5c) 
(5d) 

(5e) 



and 7(/; Qa, ^l) is the overlap reduction function, which describes the amplitude of the 
correlation of the gravitational wave signal between the two detectors as a function of their 
relative orientation ITsHT? . 



B. The overlap reduction function 

The overlap reduction function reflects frequency- dependent correlation of the gravita- 
tional wave signal in the two detectors and is expressed in terms of their relative orientation 
and separation. To express 7 for the ALLEGRO and LLO detectors, focus first on the two 
arms of the LLO detector, which define a plane (see figure |l]). Let fix be the unit vector 
along the projection onto this plane of the unit vector pointing from the LLO vertex toward 
the ALLEGRO bar's midpoint. Taking to be orthogonal to the plane and in the direc- 
tion of increasing altitude, define riy in the plane to form, with and riy a right-handed 
coordinate system. 

By construction the arms of the LLO detector also lie in this plane and, by convention, 
they are referred to as the X and Y arms. Denoting the unit vector in the direction from 
LLO vertex along the X arm as X and the unit vector in the direction from the vertex along 
the Y arm as 1^, X and Y form a right-handed coordinate system with n^. It is convenient 
to introduce the bisector of the pair {X,Y), in the direction of increasing X and Y, and 
define the angle ai to be the angle between the bisector and rix- Then, writing the TT-gauge 
gravitational wave strain as hij, the LLO detector responds to the superposition 

hi ■■= hijCtl (6) 

where 

dhic^L) ■= ^ [sin(2(TL) (n^. ® - fiy ® fiy) - cos(2(Tl) (n^ ®ny + hy® n^.)] (7) 
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of the incident gravitational waves. In computing the overlap reduction function, di char- 
acterizes the LLO detector. 

Assume that symmetry axis of the ALLEGRO bar lies in the n^. x riy plane.Q and define 

a A '■= (the angle between the ALLEGRO symmetry axis and n^,) . (8) 

In terms of cr^ the ALLEGRO detector responds to the superposition 

hA ■■= hijd'l, (9) 

where 

Aa^o-a) '■= {fix cosaA + ny sin aA) 

{fix cos (7a + fly sina a) — -I, (10) 

o 

of the incident gravitational waves. In computing the overlap reduction function, cIa char- 
acterizes the ALLEGRO detector. 

Finally, define the functions pk{a) by 




a 



—10a 5 

-lOa^ 40a -50 I | ji(a) I , (11) 




5aV2 -25a 175/2 

where the jk are the spherical Bessel functions of order k, 

a := 27r/L/c, (12) 

and L is the length of the baseline between the two detectors. The functions Pk{<y) char- 
acterize the frequency dependent part of the sensitivity of the detector pair to a stochastic 
gravitational wave background. 

Taking all these pieces together, the overlap reduction can be expressed 

7(/; Qa, ^l) ■= Pi{a)dA : cZl + p2{a) {fix ■ c^a) ■ {di ■ n-x) 

+ ps{a) {fix ■ dA ■ fix) {fix ■ di ■ fix) , (13) 

Table |I| gives the relevant geographic parameters describing LLO detectors table |T|does 
the same for the ALLEGRO detector ||14 . 



III. APPLICATION: ALLEGRO AND LLO 

The LLO detector orientation is fixed; however, the orientation of the ALLEGRO detec- 
tor may be changed by rotating ALLEGRO in its horizontal plane (cf. figure |l]). This degree 
of freedom is described by the angle a a, defined in equation R|. As a a varies, C{flA,^L) 



^This is a good approximation for the ALLEGRO/LLO detector pair. 
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will change through the dependence of 7 on cr^. To express that variation write the output 
of detector k as the sum of a gravitational wave signal hk and detector noise n^: 



Sk{t) ■.= hk{t)+nk{t). (14) 
We can then write the ensemble average of the correlation C{flA, ^l) as a function of a^: 

(15a) 



where 



{hA, hLj ■■= Tint / df 



C(o-a) := (< hA, hL> + < ua, ul >) , 
?>Hl j\f;(TA)nGw{mo,Gw{f) 



{nA,nL) :=Tint / df SalU'^cta 



207r2 /6 {SA{f)SL{f) + SAL{f?y 

7(/;(jA)fiGW,o(/) 



' P{SAU)SLU) + SALU?y 

^GW := ( The actual stochastic gravitational wave spectrum ) 



and 



SAL{f] (Ta) ■ = 



f The ALLEGRO/LLO noise cross-spectral density, 

y which in general may depend on the orientation angle a a 



(15b) 

(15c) 
(15d) 

(15e) 



Since 7 depends on the orientation a a of the ALLEGRO detector, changing ALLEGRO 's 
orientation changes C and allows us to modulate the gravitational wave contribution to C 
in a predictable way. 

Signal correlations between LLO and ALLEGRO occur only in narrow bands centered 
on the two ALLEGRO bar resonances (/<, />; cf. table |ITT|). Over the band bounded by the 
resonances the LLO noise power spectral density Sl should be approximately constant and 
we expect that ^gw will be constant as well. Additionally, for the ALLEGRO/LLO detector 
pair a (cf. Eq. |l^) is small (approximately 0.8) and does not change significantly between 
the two resonances, so that the overlap reduction function can be treated as frequency 
independent where the integrands in either of equations |15c| and |15b| are significant. Further 
assuming that Sal-, the ALLEGRO-LLO instrumental noise cross-spectral noise density, is 
independent of a a and much smaller than either Sl or Sa, we obtain 



^GW,o(/o 



/o'^a(/o)^l(/o) 



where 



A/: = 

Sa:-- 



3Hq fiGw(/o) 



207r2 



SA{f>) SA{f<) 



+ 7(/o;(^a)^al(/o) 



[SA{f>)+SA{f<)]. 



(16a) 

(16b) 
(16c) 



As a A varies the contribution of the stochastic signal to C {{Ha, hi), which is quadratic in 7) 
varies differently than the contribution of the instrumental noise {{ua, n^), which is linear in 
7). Figure shows, as a solid line, the dependence of 7 on a^. For reference, the dotted line 
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shows the dependence of 7 on a a at zero frequency. Note how 7 is approximately sinusoidal 
in 2a. Figure shows a schematic of the ALLEGRO bar orientation corresponding to the 
extrema and null of 7 as shown in figure 0. 

Since the two additive terms in equation |16a| depend on 7(/o, (Ta) differently, varying 
(Ja modulates the contribution to C of any correlated noise differently than it modulates 
the contribution of a real signal. We can use this differential modulation to eliminate the 
contribution of any correlated noise Sal that is independent of a a- Denote the angle a a for 
which 7 is maximized as crA,max; similarly, denote the angle a a for which 7 is minimized as 
c"A,min- Suppose we make an observation of duration Tint,max with ALLEGRO oriented at 
angle CTyi max, and another observation of duration Tint,min at angle a^^min- The expectation 
value of C for these two observations is 



and 



where 



A/ 



GW,0 



JqSaSl 



ml 



7max^GW + 7max'S'AL 



int,min 

^3 



GW,0 



USaSi 



7min^GW + 7min5'AL 



7niax • 7(''"A,max) 
7min • 7(''"A,min) 

The combination of these two observations 



Cn: = 



7max -^int , max ( Cmin ) 7iiiin -^int , min C* ( C^n 



7max-^int,max ~l~ |7min| -^i 



(17a) 



(17b) 



(17c) 
(17d) 



(18) 



int,min 



thus has an expectation value that is independent of the correlated noise Sal- In the present 
circumstance, 

70 • 7max — 7min (-^'^) 

If we also make the observations of equal duration, 

Tint,max -^int,min (20) 

then, following past convention and defining the signal-to-noise ratio po of the observation 
Co as the ratio of Cq to its ensemble variance we haveQ 

7o^^Gw(|/o|) 



Po 



107r2 



/oV^a(/o)5l(/o) + 5al(/o)^ 



(21) 



^Here we assume that both Sa and Sl are much greater than either \Sal\ or the corresponding 
power spectral density of the stochastic signal Sh- Were this not the case we would likely be able 
to identify the origin of the correlated noise and either isolate the detector pair from it, or regress 
it from the data during analysis. 
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IV. NUMERICAL RESULTS 



Table |IT| describes the theoretical limiting operating characteristics of the current AL- 



LEGRO [Q,0,0], the initial LIGO Livingston detector [20|, and the planned upgrades to 



ALLEGRO ||14[ and LIGO |21|. The last two rows give the theoretical limiting sensitivities 
(90% confidence bounds) for flG^j^(900Hz) that are achievable by the detector pairs LIGO I 
+ ALLEGRO and Advanced LIGO + Upgraded ALLEGRO in a one year observation. These 
limits can be reached only by identifying and accounting for non-gravitational wave inter- 
detector correlations. Unaccounted correlations of non-gravitational wave origin introduce 
a systematic error that quickly becomes the limiting factor in an upper limit determination. 

Consider, for example, the current generation of ALLEGRO and LLO detectors jointly 
observing in the presence of a correlated noise 

\SAL{fo)\ = 10-' {SASLf\ (22) 

Assume that the stochastic gravitational wave signal amplitude is much smaller: ^gw equal 
to 10~^. Suppose first that we are ignore the possibility that the correlated noise compo- 
nent (represented by Sal) may be present. Then we would leave the ALLEGRO detector 
orientation fixed in such a manner as to maximize the overlap with LLO. The dashed lines 
in figure |^A show, as a function of observing time, the 90% confidence interval (following 
the construction of [^) associated with an observed cross-correlation C (cf. eq. equal to 
the ensemble mean C. After approximately 0.25 y this most likely observation is clearly no 
longer consistent with the actual stochastic gravitational wave background amplitude, owing 
to the systematic error made by excluding the possibility of a correlated noise background. 
As the observation time increases, the confidence interval on Qqw shrinks, asymptoting on 
the amplitude of the correlated noise (Sal) interpreted as a stochastic gravitational signal. 

On the other hand, suppose we admit the possibility of a correlated noise background, of 
unknown cross-spectral density, changing the orientation of the ALLEGRO detector mid- way 



through the observation in order that we can construct Co (cf. eq. [T^), which is independent 
of Sal- Again referring to figure |^, the thin gray line shows the 90% confidence interval 
(following the construction of ^^) on Qqw when the observed Co is equal to its ensemble 
mean. The confidence interval is, in this case, always consistent with a stochastic gravita- 
tional wave background amplitude flow of 10~^. Additionally, in less than 0.45 y this 90% 
bound limits the signal amplitude to less than the correlated background noise amplitude. 

In this example the modulation technique described here provides, in approximately 
0.45 y, a bound on the stochastic signal below the correlated noise background amplitude, 
interpreted as a stochastic gravitational wave signal. Figure |^ shows the integration period 
required, using this technique, to limit the stochastic background to an amplitude less than 
the correlated noise background as a function ISALifo)]^^"^- The solid line corresponds to 
the LIGO I/ALLEGRO detector pair while the dashed line, labeled (ii), corresponds to the 
Advanced LIGO/Upgraded ALLEGRO detector pair. Since, with fixed detectors, the upper 
limit on the stochastic signal strength is always above the amplitude of the correlated noise, 
figure ^3 shows that, after 1 y of observation with LIGO I + ALLEGRO, an unaccounted 
for correlation in the background at the level of a/ SalW) ~ 3 x 10^^'^ l/\/Hz compromises 
the measurement. Similarly, after 1 y of observation with Advanced LIGO + Upgraded 
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ALLEGRO, a correlated background with a strain spectral density of ^/SalU) ~ 2 x 
10~^^ 1/V Hz will compromise a simple correlation measurement that does not account 
properly for environmental correlations. 



V. DISCUSSION 

Weber-bar gravitational wave detectors, like ALLEGRO, have placed progressively lower 
upper limits on the stochastic gravitational wave background amplitude near 900 Hz |2^-|2S[. 
The present best upper limit on ^gw; set by Astone et al. at 900 Hz p6| , p8| , used two cryo- 
genic bar detectors operating over a fairly long (600 km) baseline. During a relatively short 
measurement time they saw no excess correlation between the two bar signals. The planned 
LIGO-AUegro observations will be made over a much shorter baseline (40 Km); consequently, 
it correlated environmental noises may be more significant and greater attention will need 
to be paid to identifying and dismissing them. 

In the presence of an undetected correlated background, the best upper limit that can 
be set will eventually become limited by the bias introduced by the presence of unaccounted 
correlated detector noise. For the ALLEGRO-LLO pair, correlated noise at a level of ~ 
10~^ the geometric mean noise spectral density of the two detectors will compromise the 
observation after less than 1 year of integration time. The possibility that geophysical or 
other terrestrial correlations could compromise the observation in a single orientation of 



the detector pair are discussed in [p!q -[T8[|; however, those discussions focus on estimating 
the maximum allowable correlated background for determining a given value for an upper 
limit. We show here how rotating the ALLEGRO detector allows us to distinguish between a 
correlated background and a weak, stationary stochastic gravitational wave signal, removing 
this limit on the sensitivity of a measurement. 

The ALLEGRO detector has recently been moved to new quarters and mounted on 



an air-bearing |14], allowing its orientation with respect to LLO be modulated in manner 
described here. In a real observation campaign, the frequency with which the detector can 
be re-oriented is limited from above by the loss of observing time that is introduced through 
the disturbance of rotating of the detector and the need for the detector to settle down after 
rotation. Similarly, it is limited from below by the desire to make measurements in multiple 
orientations. A reasonable compromise would be to re-orient ALLEGRO every 3 — 5 
months, while making sure that that exact period is not commensurate with any obvious 
seasonal or annual cycles. 

The degree of improvement possible with this technique will depend on (i) whether the 
terrestrial background has any component with a characteristic quadrupolar signature that 
aliases into the stochastic gravitational wave background signature and (ii) the degree to 
which the background is stationary over the separate periods of measurement in different 
orientations. Nonetheless, an experiment that modulates the stochastic gravitational wave 
background signature will improve the quality of any long term observation of I^gw- 
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APPENDIX A: THE OPTIMAL FILTER IN THE PRESENCE OF CORRELATED 

DETECTOR NOISE 

In this appendix we evaluate the kernel Q (cf. eq. |l^) that maximizes the signal-to-noise 
ratio in equation ^, when a correlated noise component is present. 
The general two-detector cross-correlation is given by: 

rT/2 pT/2 

C{T):= dt' dt"si{t')s2{t")Q{t' -t") (Al) 

J-T/2 J-T/2 

where Si{t) are the signals from the two interferometers and Q{t) is the optimal filter 
kernel to be determined. Si{t) = hi{t) + ni{t), where hi{t) is the signal and ni{t) is the 
environmental-plus-instrumental noise in each detector. 

The correlation can also be expressed in the frequency domain. With the definition of 
the Fourier transform of a function r{t) given by: 

/oo 
dt e-^^-f' r{t), (A2) 
'OO 

the correlation measurement is also equal to, 

/oo l*CO 
df / df 5T{f - frslU)W)QU') ■ (A3) 
-oo J —oo 

Here Sj(/), and Q{f) are the Fourier transforms of the signals Sj(t) and Q(t — t'). For real 
Q{t — t'), the positive and negative frequency values of its Fourier transform are related by 
Q{—f) = Q*{f)- The function 5t(/ — /') is the finite-time approximation to the Dirac delta 
function 6{f — f) defined by 
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Hf) := r dt e-""' := 2!^. (A4) 

J-T/2 TT/ 

which arises owing to the windowing property of the finite duration (T) measurement. 
Generahzing on both the signals and the noise are now assumed to be correlated: 

ChWhif)) ■■= S{f - /')^^^^^Gw(|/|)7(l/l) (A5a) 
(^*(/)%(/')):=^5(/-/')5.,(|/|) (A5b) 

(n*(/)n,(/')<(r)%(r)) := \ {s^ms,{\fmf + rw + n+ (asc) 

s^M)s^A\f"\w-mf"-n + 

SMDSM'IW - rW - f")) 

Equation [A5c| follows from the moment theorem for real Gaussian noise p9|, Chapter 3] . 
Si{\f\) is the power spectral density of the detector noise for the i*^ detector and Sij{\f\) is 
the cross-spectral density for the pair (^, j). The factor of 1/2 in the power spectral densities 
arises from the definition of the power spectrum as a function for / > 0. 

The cross-correlated signal is a function of relative detector orientation through the 
overlap reduction factor. Call C+,_ the integrated measurement made for two orientations 
so that the contribution to each from the stochastic gravitational wave background is given 

by 

{KU)h2U)) = ±^^^^Gw(|/|)7(l/l) (A6) 

For a total measurement interval T, the strategy will be to measure first in one orientation 
for a fraction of the time, and then to measure in the other orientation for the remainder of 
the time. Consider the case in which the two intervals are each equal to T/2. The two of 
measurements have the following expectation values: 

(C^,-) ■■= T/2 df (±^^^^7(1/1) + Qif). (A7) 

The two quantities C+ _ can be used to eliminate from the expectation value of the 
derived signal the effect of the correlated detector noise: 

(C) := {C^ - C_) := T df ^^^7(1/1)) Qif) ■ (A8) 

The optimal filter function follows from maximizing the the signal-to-noise ratio of the 
cross-correlation: 

SNR := (A9) 
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where cj^ = (C^) — (C)^ is the variance of the cross-correlation signal C. Now the measure- 
ments C+^_ take place during two distinct intervals each of duration T/2. It is reasonable 
to assume that the fluctuations in these two measurements are not correlated. Then, 

4 := {C') - {Cf (AlOa) 

(AlOb) 

2cxl^ (AlOc) 

The last simplification follows from the assumed stationarity of the noise and equal mea- 
surement intervals with identical statistical properties. 

The technique used to evaluate cx^^ is similar to that of |TB|, Eq. 3.61]. First assume 
that the noise intrinsic to the two detectors is much larger in magnitude than the stochastic 
gravitational wave background. Then the variance of the measurement will be dominated 
by the detector noise and not the astrophysical signal. In this case, 

4, := {CD - {C^f (Alia) 

/oo POO poo 

df dk / dk'5T{f - f')5T{k- k') 
oo J ~oo J ~oo 

X {{nl{f)Ur)nl{k)n^{k')) - {nl{f)n^{f)) {nl{k)n^{k'))) QiDQik') (Allb) 

/OO POO POO POO 

df / df dk / dk'5Tif-nSTik-k')QinQik') 
-OO J —oo J ~oo J ~oo 

X \ [Sr{\f\)S^{\r\)5{f + k)5{r + k') + S,,{\f\)Sr,m)5{f - k')5if - k)] , (Allc) 

where Si and 5*2 are the power spectral densities of the noise in detectors 1 and 2, and Su 
is the cross-spectral density of the noise in the two detectors. 

The integrals over k, k' collapse due to the 5-functions, leaving: 

df rf/'4(/-/')(5i(l/l)52(|/'|) + 5i2(|/|)5.,(|/'|))Q(/)Q*(/0 (Al2a) 



4 
T 



oo 
oo 



df {Sims2m+Su{\f\r)\Q{f)\'. (Am) 



The last expression is obtained by approximating one of the finite-time delta functions, 
^nif — /') ; as an ordinary Dirac delta function while evaluating the other at / = /'. 
Thus the signal and its variance are given by: 

/oo / o tt2 \ 

Jf (^l/r'figw(|/|)7(l/l)jQ(/), (Al3a) 

rji POO 

4^7/ df {Si{\f\)Sim + Su{\f\nQ{f)\'- (Al3b) 

^ J — oo 

Derivation of the optimal function Q{f) proceeds along the same lines as in H. The 
analogous result is: 

nrn 7(l/l)^gw(|/|) , . 

^^^'"^ \ f\HSi{\f\)sm) + Su{\f\y)' ^ ' 
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It can be seen that the effect of correlated inter-detector noise is to replace the product of 
the single-detector power spectral densities with a sum of this product and the square of 
the contribution coming from the cross-spectrum. If the cross-spectral density is sufficiently 
small compared to the single-detector noise spectral densities, then the effect on the optimal 
filter will be correspondingly small. 
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TABLES 



TABLE L Geographic Data for LIGO Livingston Laboratory (LLO). Positions are with respect 
to the Earth Centered Frame (ECF): z pierces the earth at the north pole, x pierces the earth at 
the intersection of the prime meridian and the equator, and y = z x x 



Quantity 


Symbol 


Value 


Units 


LLO Vertex 


{Xe,Ye,Ze} 


{-74276.044, -5496283.721, 3224257.018} 


m 




{h, 4>, A} 


{-6.568, Af30° 33' 6.871", W90° 48' 50.229"} 


{m,dms,dms} 


X arm 








unit vector 


{X ■ fix, X ■ny,X ■ fiz} 


{-0.954574,-0.1415805,-0.2621887} 


ECF 


Y arm 








unit vector 


{Y ■ nx, Y fiy.Y ■ fiz} 


{0.2977412,-0.4879104,-0.8205447} 


ECF 


Bearing of 






Reference is 


ALLEGRO at 




sm.88°w 


geodetic 


LLO Vertex 






north 


Angle between 








LLO arm 






Degrees, 


bisector and 




39.59° (Bearing: 527.28° VK) 


measured 


LLO- 






CCW from 


ALLEGRO 






baseline 


baseline 









TABLE IL Geographic Data for ALLEGRO Bar Detector at LSU. For the definition of the 
Earth Centered Frame see the caption of table | 



Quantity 


Symbol 


Value 


Units 


ALLEGRO Vertex 


{Xe, Ye, Ze] 


{-113258.848, 5504077.706, 3209892.353} 


m 




{<A,A} 


{iV30° 24' 45.110", Ty91° 10' 43.766"} 


{dms, dms} 


Bearing of LLO 






Reference is 


Vertex at 




N66.67°E 


geodetic north 


ALLEGRO 








Angle between 




Correlation maximum:7max(921i/z) = 0.953 


Degrees, 


ALLEGRO bar axis 




-5.60° (Bar axis bearing: 572.08° W") 


measured 


and LLO- ALLEGRO 




Correlation null: ^nulli^'^^Hz) = 0.0 


CCW 


baseline for 




40.52° (Bar axis bearing: 526. 15°!^) 


from baseline 


various values 




Correlation minimum:7min(921//z) = —0.893 




of correlations 




84.60° (Bar axis bearing: 517.92° 




LLO - ALLEGRO 








baseline distance 


L 


42269.951 


m 


Angle subtended by 








LLO - ALLEGRO 


(3 


0.358° 


Degrees 


baseline at center 








of Earth 
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TABLE III. Operational characteristics of the LLO and ALLEGRO detectors. 



Quantity 


Lower ALLEGRO 
resonant frequency, j< 


Upper ALLEGRO 
resonant frequency, /> 


Units 


Frequency 


896.8 


920.3 


Hz 


ALLEGRO 
bandwidth, Af 


1 


1 


Hz 


Upgraded ALLEGRO 
bandwidth |14], A/ 


50 


50 


Hz 


ALLEGRO sensitivity, 

Hf) 


1.8 X 10-21 


0.85 X 10-21 


1/VHz 


LIGO I sensitivity |2|, 

Hf) 


1 X 10-22 


1 X 10-22 


1/VHz 


Adv. LIGO narrowband 
spnsitivitv II2TII h( f) 


2 X 10-24 


2 X 10-24 


1/VHz 


l^min for LIGO I + 

ALLEGRO after 1 year 


1 X IQ-i 


at 90% confidence 


^min for adv. LIGO + 
Upgraded ALLEGRO 
after 1 year 


2.9 X 10-4 


at 90% confidence 



17 



FIGURES 




FIG. 1. A schematic diagram showing how we characterize ALLEGRO and LLO orientations 
with respect to geodetic north and the LLO-to- ALLEGRO baseline. 
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f = Hz 

/ 



Geodetic 


North 




LLO/ALLEGRO 




Rparinn 




^ S66.68 W 


Tmax: S72.08W J 




Ynul|:S26.15W 


ymin:S17.92E 




FIG. 2. The overlap reduction function 7, which characterizes the dependence of the ALLE- 
GRO-LLO correlation function to an isotropic stochastic signal, depends on the relative orientation 
of the two detectors. Here we show how this function varies with cja, the angle between the AL- 
LEGRO bar axis and the LLO/ALLEGRO baseline (cf. figure |l|). The bold, solid line shows the 
variation of 7 with a a at the operating frequency of the ALLEGRO detector; for comparison, 
the dashed line shows the same quantity at DC. Inset B shows the orientation of the ALLEGRO 
relative to the LLO/ALLEGRO baseline, when 7 vanishes (7nuii) and takes on its minimum (7min) 
and maximum (7max) values. 
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FIG. 3. A) The dashed lines mark the expected 90% confidence interval, as a function of 
the observing time, on a stochastic gravitational wave background when a much larger, but unac- 
counted for, correlated terrestrial noise source is present with a cross spectral density amplitude 
just 10~^ the (geometric) mean noise power spectral density in the ALLEGRO and LIGO I detec- 
tors. The heavy solid line is the amplitude of the terrestrial noise, (mis) interpreted as a stochastic 
gravitational wave signal. Note how, after approximately 3 months, the observations are no longer 
consistent with stochastic gravitational wave amplitude significantly less than the amplitude of the 
correlated terrestrial noise. The thin line marks the upper limit on the stochastic signal, again as 
a function of time, when the modulation technique described in this paper is used to make the 
measurement. The measurement is no longer biased by the terrestrial noise and the upper limit is 
less than the correlated terrestrial noise amplitude, in this example, in 0.45 y. B) The integration 
time needed for the upper limit, estimated by the modulation technique described here, to be less 
than the amplitude of the correlated terrestrial noise amplitude (i.e., to reach the crossing point 
marked by the bold circle in panel A) as a function of the cross-spectral density. Curve (i) is for 
LIGO I + ALLEGRO; (ii) is for advanced LIGO + Upgraded ALLEGRO (cf. table |lT|). 
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